Quantum Phase Transition of a Magnet in a Spin Bath 
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The excitation spectrum of a model magnetic system, LiHoF 4 , has been studied 
using neutron spectroscopy as the system is tuned to its quantum critical point by 
an applied magnetic field. The electronic mode softening expected for a quantum 
phase transition is forestalled by hyperfine coupling to the nuclear spins. We show 
that interactions with the nuclear spin bath control the length scale over which the 
excitations can be entangled. This generic result limits how far it is possible to 
approach intrinsic electronic quantum criticality. 

The preparation and preservation of entangled quantum states is particularly relevant for the 
development of quantum computers, where interacting qubits must produce states sufficiently 
long-lived for meaningful manipulation. The state lifetime, typically referred to as a decoher- 
ence time, is derived from coupling to the background environment. For solid state quantum 
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computing schemes, the qubits are typically electron spins, and they couple to two generic 
background environments 0. The oscillator-bath — i.e. delocalised environmental modes (0), 
such as thermal vibrations coupled via magneto-elastic terms to the spins — can be escaped by 
lowering the temperature to a point where the lattice is essentially frozen. Coupling to local de- 
grees of freedom, such as nuclear magnetic moments which form a spin-bath, may prove more 
difficult to avoid, as all spin-based candidate materials for quantum computation have at least 
one naturally occurring isotope that carries nuclear spin. Experimental work on this subject has 
been largely restricted to the relaxation of single, weakly interacting, magnetic moments such 
as those on large molecules © while much less is known about spins as they might interact in a 
real quantum computer. In this regard, the insight that quantum phase transitions are a good 
arena for looking at fundamental quantum properties of strongly interacting spins turns out to 
be valuable, as it has already been for explorations of entanglement. In particular, we show 
that coupling to a nuclear spin bath limits the distance over which quantum mechanical mixing 
affects the electron spin dynamics. 

Quantum phase transitions (QPTs) are transitions between different ground states driven 
not by thermal fluctuations, but by quantum fluctuations controlled by some parameter such as 
doping, pressure and magnetic field (j5]|6|). Much of the interest in QPTs stems from their im- 
portance for understanding materials with unconventional properties - heavy Fermion systems 
and high temperature superconductors. However, these materials are rather complex, and do 
not easily lend themselves to a universal understanding of QPTs. To this end it is desirable 
to identify quantum critical systems with a well-defined and solvable Hamiltonian, and with a 
precisely controllable tuning parameter. One very simple model displaying a QPT is the Ising 
ferromagnet in a transverse magnetic field ([7||5]|P|E1>, with the Hamiltonian: 

;>t (i) 

ij i 

In the absence of a magnetic field, the ensemble of the two degenerate o z = ±1 states orders 
ferromagnetically below a critical temperature T c . The transverse field Y mixes the two states 
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and even at zero temperature leads to destruction of long-range order in a QPT at a critical 
field r c . In the ferromagnetic state at zero field and temperature, the excitation spectrum is 
momentum independent and centered at the energy 4 J2j 3ij associated with single spin reversal. 
Upon application of a magnetic field, however, the excitations acquire a dispersion, softening 
to zero at the zone centre q = when the QPT is reached. 

We investigate the excitation spectrum around the QPT in LiHoF 4 , which is an excellent 
physical realisation of the transverse field Ising model, with an added term accounting for the 
hyperfine coupling between electronic and nuclear moments (|771l70l U2j) . The dilution series 
LiHo. r Y 1 _. I .F 4 is the host for a wide variety of collective quantum effects, ranging from tun- 
nelling of single moments and domain walls to quantum annealing, entanglement and Rabi 
oscillations il3\\14l\T5\\16\\17b . These intriguing properties rely largely on the ability of a trans- 
verse field, whether applied externally or generated internally by the off-diagonal part of the 
magnetic dipolar interaction, to mix two degenerate crystal field states of each Ho ion. The Ho 
ions in LiHoF 4 are placed on a tetragonal Scheelite lattice with parameters a = 5.175 A and 
c = 10.75 A. The crystal-field ground state is a T 3 4 doublet with only a c component to the 
angular momentum and hence can be represented by the o z = ±1 Ising states. A transverse 
field in the a-b plane mixes the higher-lying states with the ground state producing a splitting 
of the doublet, equivalent to an effective Ising model field. The phase diagram of LiHoF 4 (Fig. 
[l]\) was determined earlier by susceptibility measurements UTk. and displays a zero field T c of 
1.53 K, and a critical field of H c = 49.5 kOe in the zero temperature limit. The same measure- 
ments confirmed the strong Ising anisotropy, with longitudinal and transverse ^-factors differing 
by a ratio of 18 (I7DI) . The sudden increase in H c below 400 mK was explained by alignment of 
the Ho nuclear moments through the hyperfine coupling. Corrections to phase diagrams due to 
hyperfine couplings have a long history (ITffl) . and were noted for the Lii?£F 4 (RE=Rare earth) 
series, of which LiHoF 4 is a member, over twenty years ago (19). What is new here is that ap- 
plying a transverse field and employing high resolution neutron scattering spectroscopy allows 
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us to carefully study the dynamics as we tune through the quantum critical point. 

We measured the magnetic excitation spectrum of LiHoF 4 using the TAS7 neutron spec- 
trometer at Ris0 National Laboratory, with an energy resolution (full width at half maximum) of 
0.06-0. 1 8 meV SWb . The transverse field was aligned to better than 0.35° and the sample cooled 
in a dilution refrigerator. At the base temperature of 0.31 K, giving a critical field of 42.4 kOe, 
the excitation spectrum was mapped out below, at and above the critical field (Fig. For 
all fields a single excitation branch disperses upwards from a minimum gap at (2,0,0) towards 
(1,0,0). From (1,0,0) to (1,0,1), the mode shows little dispersion but appears to broaden. The 
discontinuity on approaching respectively (1, 0, 1— e) and (1+e, 0, 1) reflects the anisotropy and 
long-range nature of the magnetic dipole coupling. However, the most important observation 
is that the (2,0,0) energy, which is always lower than the calculated single ion energy (~ 0.39 
meV at 42.4 kOe), shrinks on increasing the field from 36 to 42.4 kOe, and then hardens again 
at 60 kOe. At this qualitative level, what we see agrees with the mode softening predicted for 
the simple Ising model in a transverse field. However, it appears that the mode softening is 
incomplete. At the critical field of 42.4 kOe it retains a finite energy of 0.24 ± 0.01 meV. Fig. 
[TjB, showing the gap energy as a function of the external field, makes this result very apparent. 

To obtain a quantitative understanding of our experiments we consider the full rare-earth 
Hamiltonian, which closely resembles that of HoF 3 (|271l2"2|) . Each Ho-ion is subject to the 
crystal field, the Zeeman and the hyperfine coupling. The interaction between moments is dom- 
inated by the long-range dipole coupling, with a small nearest neighbour exchange interaction 

i 

n.n. 

ij af3 ij 

where J and I are respectively the electronic and nuclear moments, and for 165 Ho 3+ J = 8 
and I = 7/2. Hyperfine resonance (23) and heat capacity measurements (T2~?t show that the 
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hyperfine coupling parameter A = 3.36 yueV as for the isolated ion, with negligible nuclear- 
quadrupole coupling. The Zeeman term is reduced by the demagnetization field. The normal- 
ized dipole tensor D a/3 (ij) is directly calculable, and the dipole coupling strength is simply fixed 
by lattice constants and the magnetic moments of the ions at 3d = (9^b) 2 N = 1.1654 /ieV. 
This leaves as free parameters various numbers appearing in the crystal field Hamiltonian JicF 
and the exchange constant # 12 . The former are determined d25t largely from electron spin res- 
onance for dilute Ho atoms substituted for Y in LiYF 4 , while the latter is constrained by the 
phase diagram determined earlier (TTOt (Fig.l A). We have used an effective medium theory © 
previously applied to HoF 3 (26) to fit the phase diagram, and conclude that a good overall de- 
scription, except for a modest (14%) overestimate of the zero field transition temperature is 
obtained for # 12 = — 0.1 yueV. Based on quantum Monte Carlo simulation data others d271) have 
also concluded that $ 12 is substantially smaller than do- 

Having established a good parametrization of the Hamiltonian, we model the dynamics, 
where expansion to order 1/z (z is the number of nearest neighbors of an ion in the lattice) 
leads to an energy dependent renormalization [1 + £(a;)] -1 (of the order 10%) of the dynamic 
susceptibility calculated in the random phase approximation (RPA), with S(o;) evaluated as 
described in (l2oT) . For the three fields investigated in detail, the dispersion measured by neutron 
scattering is closely reproduced throughout the Brillouin zone. As indicated by the solid lines in 
Fig. |3 the agreement becomes excellent if the calculated excitation energies are multiplied by 
a renormalization factor Z = 1.15. The point is not that the calculation is imperfect, but rather 
that it matches the data as closely as it does. Indeed, it also predicts a weak mode splitting of 
about 0.08 meV at (1, 0, 1 — e), consistent with the increased width in the measurements. The 
agreement for the discontinuous jump between (1,0, 1 — e) and (l + e,0, 1) due to the long-range 
nature of the dipole coupling, shows that this is indeed the dominant coupling. 

Fig. [l)B-C illustrates the simple origin of the incomplete softening and enhanced critical 
field, which is easiest to understand if we start from the polarised paramagnetic state above H c , 
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where the experiment, the purely electronic calculation, and the theory including the hyperfine 
coupling all coincide. At high fields the only effect of the hyperfine term is to split both the 
ground state and the electronic excitation modes into multiplets which are simply the direct 
products of the electronic and nuclear levels, with a total span of 2A(J)I ~ 0.1 meV as illus- 
trated in Fig.[l]C. Upon lowering the field, the electronic mode softens and would in the absence 
of hyperfine coupling reach zero energy at H® = 36 kOe. The hyperfine coupling, however, 
mixes the original ground and excited (soft mode) states already above H c . As this happens, 
the formation of a composite spin from mixed nuclear and electronic contributions immediately 
stabilizes ordering along the c-axis of the crystal. In other words, the hyperfine coupling shunts 
the electronic mode, raising the critical field to the observed H C =42A kOe, where the mode 
reaches a non-zero minimum. This process is accompanied by transfer of intensity from the 
magnetic excitation of electronic origin to much lower energy (in the 10 fieV range) soft modes 
of entangled nuclear/electronic character. Cooling to very low temperatures would reveal these 
modes as propagating and softening to zero at the quantum critical point, but at the temperatures 
reachable in our measurements there is thermalization, dephasing the composite modes to yield 
the strong quasi-elastic scattering appearing around Q = (2, 0, 0) and zero energy at the critical 
field, as in Fig. |2j 

The intensities of the excitations are simply proportional to the matrix elements | (/ 1 Yjj exp(iQ- 
Rj)Jj~\0)\ 2 , and therefore provide a direct measure of the wavefunctions via the interference 
effects implicit in the spatial Fourier transform of Fig. |3] shows intensities recorded along 
(h, 0, 0) for the three fields 36, 42.4, and 60 kOe. They follow a momentum dependence char- 
acterized by a broad peak near (2,0,0), which is well described by our theory. In the absence of 
hyperfine interactions, the intensity at H® would diverge as q approaches (2, 0, 0), reflecting that 
the real-space dynamical coherence length £ c of the excited state grows to infinity. The finite 
width of the peak observed at H c , corresponds in real space to a distance of order the inter- 
holmium spacing, and implies that because they forestall the softening of the electronic mode, 



6 



the hyperfine interactions also limit the distance over which the electronic wavefunctions can 
be entangled (4). Thus, Fig.|3]is a direct demonstration of the limitation of quantum coherence 
in space via coupling to a nuclear spin bath. The dynamical length £ c is obtained from a sum 
over matrix elements connecting the ground state to a particular set of excited states, while the 
thermodynamic correlation length £ t is derived from the equal time correlation function S(r) 
which is the sum over all final states. £ 4 diverges at second order transitions such as those in 
L1H0F4, where the quasielastic component seen in our data dominates the long-distance be- 
haviour of S(r) at T C (H). It is the electronic mode and hence £ c that dictates to which extent 
LiHoF 4 can be characterised and potentially exploited as a realisation of the ideal transverse 
field Ising model. 

Beyond providing a quantitative understanding of the excitations near the quantum critical 
point of a model experimental system, we obtain new insight by bringing together the older 
knowledge from rare-earth magnetism and the contemporary ideas of entanglement, qubits and 
decoherence. While the notion of the spin-bath was developed to address decoherence in lo- 
calised magnetic clusters and molecules (1 ), our work discloses its significance for quantum 
phase transitions. Notably, we establish that the spin bath is a generic feature that will limit 
how far it is possible to observe intrinsic electronic quantum criticality. While this may not be 
a significant factor when thinking about transition metal oxides with very large exchange con- 
stants, it could matter for rare earth and actinide intermetallic compounds which show currently 
unexplained cross-overs to novel behaviours at low (< 1 K) temperatures (see e.g. fiZBk). 

For magnetic clusters, decoherence can be minimised in a window between the oscillator- 
bath dominated high temperature and the spin-bath dominated low-temperature regions 
Our calculations suggest that the dense quantum critical magnet shows analogous behaviour. 
Here the interacting electron spins themselves constitute the oscillator bath, and the extent to 
which the magnetic excitation softens at T C (H), measured by the ratio of the zone center energy 
E c to the field-induced single ion splitting A (Fig. [TJ)), gauges the electronic decoherence. 
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E c /A achieves its minimum not at T = 0, but rather at an intermediate temperature T ~ 1 K, 
exactly where the phase boundary in Fig. UK begins to be affected by the nuclear hyperfine 
interactions. 
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Figure 1 : A Phase diagram of LiHoF 4 as a function of transverse magnetic field and temperature 
from susceptibility U(k (circles) and neutron scattering (squares) measurements. Lines are 1/z 
calculations with (solid) and without (dashed) hyperfine interaction. Horizontal dashed guide 
marks the temperature 0.31 K at which inelastic neutron measurements were performed. B 
Field dependence of the lowest excitation energy in LiHoF 4 measured at Q = (1 + e, 0, 1). 
Lines are calculated energies scaled by Z = 1.15 respectively with (solid) and without (dashed) 
hyperfine coupling. The dashed vertical guides show how in either case the minimum energy 
occurs at the field of the transition (c.f. A). C Schematic of electronic (blue) and nuclear (red) 
levels as the transverse field is lowered towards the QCP. Neglecting the nuclear spins, the 
electronic transition (light blue arrow) would soften all the way to zero energy. Hyperfine 
coupling creates a non-degenerate multiplet around each electronic state. The QCP now occurs 
when the excited state multiplet through level repulsion squeezes the collective mode of the 
ground state multiplet to zero energy, hence forestalling complete softening of the electronic 
mode. Of course, the true ground and excited states are collective modes of many Ho ions and 
should be classified in momentum space. D Calculated ratio E c /A of the minimum excitation 
energy E c to the single ion splitting A at the critical field as a function of temperature. This 
measures how far the electronic system is from the coherent limit, for which E c / A = 0. 
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Figure 2: Pseudocolor representation of the inelastic neutron scattering intensity for LiHoF 4 
at T = 0.31 K observed along the reciprocal space trace (2, 0, 0) — > (1, 0, 0) — > (1, 0, 1) — > 
(1.15,0, 1). White lines show the 1/z calculation for the excitation energies as described in 
the text. White ellipses around the (2,0,0) Bragg peak indicate 5 times the full-width-half-max 
resolution tail. 
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Figure 3: Measured intensities of the excitations along Q = (h, 0, 0) at the same values of the 
field as in Fig. 2. Lines are calculated with geometric and resolution corrections applied to 
allow comparison to the neutron data. 
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